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The coupled interaction between an ablating surface and the surrounding aerothermal 
environment is studied. An equilibrium ablation model is coupled to the LAURA flowfield 
solver, which allows the char ablation rate (m c ) to be computed as part of the flowfield 
solution. The wall temperature ( T w ) and pyrolysis ablation rate (rh g ) may be specified by 
the user, obtained from the steady-state ablation approximation, or computed from a a 
material response code. A 32 species thermo chemical nonequilibrium flowfield model is 
applied, which permits the treatment of C, H, O, N, and Si containing species. Coupled 
ablation cases relevant to the Orion heatshield are studied. These consist of diffusion- 
limited oxidation cases with Avcoat as the ablation material. The m c values predicted 
from the developed coupled ablation analysis were compared with those obtained from a 
typical uncoupled ablation analysis. The coupled results were found to be as much as 50% 
greater than the uncoupled values. This is shown to be a result of the cumulative effect of 
the two fundamental approximations inherent in the uncoupled analysis. 

I. Nomenclature 

A Blowing reduction constant = 0.5 

6 Char emissivity 

p c Char density (kg/m 2 ) 

p v Density of virgin material (kg/m 2 ) 

a Stefan-Boltzmann constant = 5.66961e-8 W/(m 2 K 4 ) 

r Shear stress 

A + Mixing length damping factor 

B 0 Blowing parameter = tti/Ch q 

B' c Non-dimensional char ablation rate = m c / C h 

B' g Non-dimensional pyrolysis ablation rate = ihg/Cn 

c c? fc Mass fraction of element k in the char gas 

c e , fe Mass fraction of element k at the boundary layer edge 

c g ,k Mass fraction of element k in the pyrolysis gas 

c w k Mass fraction of element k at the wall 

Cw,c Mass fraction of atomic carbon at the wall 

Ch Dimensional heat- transfer coefficient (kg/(m 2 s)) 

Ht Flowfield total enthalpy ( J /kg) 

h w Enthalpy of gas at wall (J/kg) 

h c Enthalpy of solid char at T w ( J /kg) 

h g Enthalpy of pyrolysis gas in equilibrium at T w and p w (J/kg) 

J k Diffusion mass flux of element k (kg/ (m 2 s) 

K c c Equilibrium constant for the heterogenous reaction C (solid) — > C(gas) (kg-mol/m 3 ) 

ra c Mass flux of surface char (kg/ (m 2 s)) 
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rn g Mass flux of pyrolysis gas at the surface (kg/(m 2 s)) 

rn Total mass flux at the surface = m c + rn g (kg/(m 2 s)) 

Me Molecular weight of carbon = 12.01 kg/kg-mol 

q c Convective heat flux at surface(W/cm 2 ) 

q c o Non-ablating convective heat flux at surface 

T w Wall temperature (K) 

Normal velocity at wall divided by friction velocity 


II. Introduction 

The assessment and design of reentry heat shields typically involve an uncoupled ablation-flowfield anal- 
ysis. This uncoupled analysis consists of applying a non- ablating flowfield prediction, which defines the 
heat-transfer coefficient, wall enthalpy and wall pressure, to an equilibrium ablation model, which computes 
the char- ablation rate (ra c ), pyrolysis ablation-rate (m 5 ), wall- temperature ( T w ), and in-depth material 
properties. The analysis is referred to as “uncoupled” because the influence of ra c , m g , and T w on the 
flowfield prediction are treated approximately within the ablation model, and hence the ablation model is 
not coupled with the flowfield model. For the diffusion-limited oxidation regime expected for Orion, the 
fundamental influence of coupling involves the prediction of m c . The prediction of rh c is sensitive to the 
chemical composition at the surface, and vice versa. This coupling requires the detailed treatment of convec- 
tion and diffusion of chemical species towards and away from the wall. For the uncoupled ablation analysis, 
two fundamental approximations are required to model this effect, assuming equilibrium ablation. These are 
defined as follows: 

1. Approximation #1 : 

The influence of m c and rn g on the heat-transfer coefficient is approximated as 

2AB 0 


Ch = Ch, o- 


(i) 


exp(2XBo) — 1 

where Ch , o is the non-abating heat-transfer coefficient. Recall that Ch is related to q c as 

q c = Ch(Ht — h w ) (2) 

This approximation has been assessed by Thompson and Gnoffo 1 and Martinelli et al. 2 for perfect gas 
flows. 

2. Approximation #2\ 

The elemental diffusion mass flux at the surface is written as 3 


Jk — (Cw,k Ce,k) (3) 

where it is assumed that Cm = Ch • This approximation allows the elemental mass balance at the wall 
to be solved algebraically for the elemental mass fractions at the wall: 3 

_ “ 1 “ B c C c ,k B g C g ,k ( .. 

= 1 + B' + B' (4) 

A discussion of this approximation is presented by Bianchi et al. 4 for graphite ablation. 

The purpose of this paper is to assess the validity of the uncoupled ablation analysis, through comparisons 
with a coupled analysis, for conditions relevant to NASA’s Orion CEV. The coupled analysis developed in 
this work is discussed in Section III, and compared with results of previous studies in Section IV. An 
assessment of Approximations #1 and #2 is presented in Section V assuming fixed values for ra c , m^, and 
T w , wdiich allows for the influence of the approximations to be clearly interpreted. Section VI compares the 
m c computed with an uncoupled and coupled analysis, with fixed values for rh g and T w . This comparison 
addresses the fundamental differences between a coupled and uncoupled analysis for the diffusion-limited 
oxidation regime of interest for Orion. 
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III. Coupled Ablation in LAURA 


This section presents the details of the developed coupled ablation analysis. The basic flowfield para- 
maters are presented in Part A, while the governing surface equations are presented in Part B. Definitions 
of the three developed solution approaches are provided in Part C. Details of the solution procedure are also 
discussed. Finally, Parts D and E discuss the applied diffusion model and turbulence model, respectively. 


A. Applied Flowfield Parameters 

For the cases presented in this work, a two- temperature thermochemcial nonequilibrium flowfield model 
was applied. The following 32 species were treated throughout the flowfield (including the wall): N, N+, 
NO, NO+, N 2 , N 2 +, O, 0+ 0 2 , 0 2 +, e, C, C+, CO, C0 2 , C 2 , C 3 , C 4 , C 5 , C 2 H, C 2 H 2 , CN, H, H 2 , 
HCN, CH, Si, SiO, Si0 2 , SiC, Si 2 , SiN. The thermodynamic properties for these species were obtained from 
Gordon and Mcbride 5 . The transport properties were obtained from Wright et al. 6,7 where available. The 
remaining species were treated using the approximate approach of Svehla 8 modified as suggested by Park. 9 
The chemical reaction rates were compiled mostly from previous Earth, 9 Mars, 10 Titan, 11 and Venus 12 
studies, with the exception of the Si rates, which were taken from the work of Mick et al. 13,14 If a given 
rate (with a different value) was available from more than one source, the hierarchy for choosing a rate was 
as follows: Fujita et al., 12 Gocken, 11 Park et al., 9 and Park. 10 Fujita et al. performed an up-to-date critical 
analysis of many of the rates, which is the reason their rate model is given preference over the others. 


B. Governing Surface Equations 

The coupled ablation capability in LAURA, assuming specified values for m c , m g , and T w , was presented 
by Gnoffo, et al. 15 . The present work extends this capability by allowing m c to be computed as part of the 
flowfield solution. The additional equation required for m c is the equilibrium char constraint: 


Cw,C 


K ( 


c,C 


(5) 


which assumes that the solid carbon char is in chemical equilibrium with the gas at the surface. The solution 
of Eq. (5) requires the atomic carbon mass fraction at the wall (c w ,c), which is assumed to be in chemical 
equilibrium at the wall temperature, pressure, and elemental composition. The elemental composition is 
obtained from the wall elemental mass balance equation: 


^hc(Cc,fc Cw,k) tilg(Cg,k Cw,k) Jk — 0 

The wall pressure ( p w ) is obtained from the normal momentum equation: 

dj)yj dv w 

; “b Pw^w 7 = t) 

dz dz 

where the normal velocity ( v w ) is obtained from the mass continuity equation: 


(6) 

(7) 


pwVw = rhc + rrig (8) 

If the pyrolysis ablation rate (rh g ) and wall temperature ( T w ) are specified, then Eqs. (5) - (8) may be 
combined with a chemical equilibrium solver to obtain ra c . 

For the present study, rh g and T w are either specified or obtained using the steady-state ablation approx- 
imation. The steady-state ablation approximation requires the solution of an approximate surface energy 
equation: 

Qc QQrad + €&T w -j- (jTl c -f- TTl g ^h w — 0 (9) 

which assumes that the char surface and virgin material recede at the same linear rate. The steady-state 
ablation approximation provides the following relationship for m g : 

rhg = (— — 1 )m c (10) 

Pc 
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C. Definitions of Applied Approaches 


To make clear the influence of Approximations #1 and #2, solutions will be presented that include both 
approximations, only Approximation #1, and neither approximation. These three different approaches are 
defined as follows: 

1. Uncoupled: This approach applies both approximations, and is identical to the typical “uncoupled” 
approach described in the Introduction. It is applied as a post-processing step to a non-ablating 
solution, obtained assuming a super-catalytic, radiative equilibrium wall. The non- ablating solution 
provides q c , o, h w , and p w . For a given ra c , m g , and T w , the elemental composition at the wall is 
obtained from Eq. (4). With this elemental composition, the species mass fractions and enthalpy 
are obtained using LAURA’s equilibrium chemistry module at the given T w and p w . This process is 
repeated while iterating on m c (and T w if it is not specified) until Eq. (5) (and Eq. (9) for T w ) are 
satisfied. Note that this approach is identical to the application of B’ tables often applied in material 
response codes, although the present approach does not actually apply any pre-computed tables. 

2. Partially- Coupled: This approach removes Approximation #1, but contains Approximation #2. It 
consists of a flowfield with coupled ablation, meaning that the injection of ablation products is treated 
in the flowfield solution. Therefore, Approximation #1 is not applied because the influence of m on 
Cjj is explicitly treated. However, this approach uses the approximation for Jk in Eq. (3), which 
means that Approximation #2 is applied. The species mass fractions and enthalpy at the wall are 
obtained identically to the uncoupled approach, except that Eq. (1) is not applied for Ch • Instead, 
Ch is computed directly from the q c computed in the flowfield solution. The iteration procedure for 
rh c and T w is also identical to the uncoupled approach. 

The application of this approach typically begins by first obtaining a non-ablating flowfield solution. 
From this, the procedure outlined in the previous paragraph to obtain ra c , T w , and c w ^i is executed 
every 5000 flowfield iterations. In between these surface computations, the pseudo-cells at the wall are 
updated every 50 flowfield iterations to maintain the computed surface properties, which are assumed 
equal to the average of the pseudo-cell values and the values in the first cell above the wall. This 
procedure is very robust and converges for a wide range of ablation rates. 

3. Fully- Coupled: This approach avoids the use of Approximations #1 and #2. It therefore provides 
the most rigorous possible application of equilibrium ablation. For this approach, Jk, w is computed in 
detail using the concentration gradients and diffusion coefficients at the wall, instead of using Eq. (3) 
as done in the Partially-Coupled approach. As a result, Eq. (6) may not be solved algebraically for 
c w k to obtain an equation similar to Eq. (4). Instead, Eq. (6) is solved using the approach developed 
by Gnoffo et ah, 15 which couples the chemical equilibrium equations with the species conservation 
equations at the surface. This approach was extended in the present study by adding Eq. (5) to this 
set of equations solved at the surface, which allows m c to be obtained as part of the surface iteration 
procedure. 

The application of this approach typically begins by first obtaining a Partialy- Coupled solution. It was 
found that the coupled set of surface equations require an reasonable initial solution for convergence. 
If the solution of T w is also desired (meaning it is not specified), then the surface energy equation is 
solved every 1000 flowfield iterations. It was found that the surface energy equation causes serious 
convergence problems if it is coupled directly to the other surface equations. 

D. Diffusion Modeling 

The “approximate-corrected” diffusion model proposed by Sutton and Gnoffo 16 is typically applied for 
aerothermo dynamic simulations. This model has been shown to produce q c values that compare very well 
with the exact model provided by the Stefan- Maxwell equations. 17 ’ 18 The species compositions through 
the boundary layer predicted by the approximate-corrected approach, however, may be noticeably different 
than the Stefan-Maxwell result. To illustrate this, Fig. 1 compares the mass fractions of elemental oxygen 
along the stagnation-line for a 3.6 m radius sphere in air with free-stream velocity and density values of 8.0 
km/s and 3xl0 -3 kg/m 3 , respectively. Note that nitrogen is the only other element present in this case, 
and its mass fraction is equal to 1 - co- The wall is assumed to be equilibrium catalytic. It is seen that the 
Stefan-Maxwell approach predicts 8% more elemental oxygen at the surface than the approximate-corrected 


4 of 15 


American Institute of Aeronautics and Astronautics 



approach. Both approaches predict a stagnation point q c of 99.9 W/cm 2 . The deviation of the elemental 
oxygen mass fractions in the boundary layer from the free-stream (and boundary layer edge) value of 0.24 
has been observed by other researchers. 19,20 

With the introduction of ablation, the deviation in the elemental mass fraction predicted by the approximate- 
corrected and Stefan-Maxwell approaches decreases. This is shown in Fig. 2, which compares all the elements 
present for an Avcoat ablation case, except nitrogen. The details of the ablator composition and ablation 
rates for this case will be present later in Fig. 13. The free-stream conditions are identical to those for Fig. 1. 
The presence of ablation tends to decrease the severe species gradients present in a pure air boundary layer. 
This is a result of nearly all the elemental oxygen combining with C to produce CO, which occurs rather 
gradually throughout the boundary layer. With all the elemental oxygen going to CO, recombination of O 2 
near the wall (and its associated large O 2 gradient) is prevented. 

The ultimate goal of the present work is to compute rh c from Eq. (5) as part of the flowfield solution. 
Since Eq. (5) is dependent on the atomic carbon mass fraction at the wall, the prediction of ra c will be very 
sensitive to the species diffusion near the surface. Although Fig. 2 showed that the approximate-corrected 
approach agrees well with the Stefan-Maxwell approach, it was found to under- predict ra c by up to 10% 
(this applies to the Fully-Coupled ablation approach). Therefore, it was decided to apply the Stefan-Maxwell 
approach for the coupled ablation cases presented in this paper. 



Figure 1: Elemental mass fractions of oxygen through the stagnation 
line boundary layer 

E. Turbulence Modeling 

Turbulence is treated in the present study using the 
Cebeci- Smith model, modified for a pressure gradi- 
ent and surface mass transfer. These modifications 
are implemented through the mixing length damp- 
ing factor, A 4 *, which is written presently as: 

A + = 26(|r w /r|) 1/2 exp(— 5.9t;+) (11) 

The term (It^/tI) 1 / 2 accounts for a pressure gradi- 
ent,, 21 and is standard in LAURA’s Cebeci-Smith 
and Baldwin- Lomax models. 22 The term exp(- 
5.9u+) accounts for surface mass transfer, as pro- 
posed by Cebeci 23 and based on the work of 
Bushnell and Beckwith. 24 The turbulent Schmidt 
(Scturb) and Prandtl ( Prturb ) numbers are set equal 
to 1.0. To show the influence of the two modifica- 
tions in Eq. (8), both a non- ablating and ablating 
case were studied with and without these modifica- 
tions. The geometry and free stream conditions are 



Figure 2: Elemental mass fractions of oxygen through the stagnation 
line boundary layer 



2 3 

s(m) 


Figure 3: Specified ablation rates and wall temperature. 
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the same as those applied in Section III.D. For the ablating case, m g , ra c , and T w are specified to the values 
shown in Fig. 3. The resulting convective heating distributions are shown in Fig. 4 for the ablating case 
and Fig. 5 for the non-ablating case. The “Baseline” result shown in these figures refers to the applica- 
tion of Eq. (11), which contains both the pressure gradient and mass-transfer modifications. The “A+ = 
26” result removes both of these modifications, while the “No mass-transfer corr.” result removes only the 
mass-transfer correction. For the ablating case, Fig. 4 indicates that removing the mass-transfer correction 
provides a slight decrease in q c , while the pressure gradient correction has a negligible influence. For the 
non-ablaitng case, Fig. 5 shows that removing the pressure gradient correction (represented by the A += 26 
result because the mass-transfer correction has no influence) provides a slight increase in q c . 




Figure 4: Convective heating along the body for ablation rates and Figure 5: Convective heating along the body for the non-ablating 

T u , specified in Fig. 3. case. 

The influence of Sc tU rb is shown in Figs. 4 and 5 through the “Turb. Sc# = 0.5” result, which applies 
Eq. (11) with a Sct^ = 0.5, instead of the value of 1.0 applied for the other results. It is clear from the 
figures that lowering the Sc tU rb significantly increases q c in the turbulent regions. Note that q c at s = 2.0 m 
is slightly larger for the ablating case with Sc tU rb = 0.5. This is a result of the increased contribution from 
the diffusive component of q c , which actually increases with the introduction of ablation. This behavior will 
be discussed in more detail later. 

The result of the Baldwin-Lomax turbulence model is also shown in Figs. 4 and 5. This model includes 
the pressure gradient correction, but not the mass-transfer correction. It is seen to compare relatively well 
with the “Baseline” model for both the ablating and non-abalting case. 

IV. Comparison with Previous Studies 

This section compares the developed coupled ablation capability with the results of three previous stud- 
ies. 25-27 These studies were chosen because they consider cases in the diffusion- limited oxidation regime of 
present interest, and they assume equilibrium ablation. As will be shown, the agreement obtained with these 
studies is not generally good. The possible reasons for this are discussed. 

A. Comparison with Results of Bartlett et al. for Avcoat 

Bartlett et al. 25 present coupled steady-state ablation results (meaning T w and rh g 
are computed from Eqs. (9) and (10)) for the stagnation point of a 3.96 m radius 
sphere. These results were obtained using the BLIMP program. 28 The ablation 
material considered was Avcoat. Both coking and mechanical silica removal were 
assumed, which leads to the elemental mass fractions shown in Table 1. 

The present Fully-Coupled results are compared with the coupled results of 
Bartlett et al. in Table 2 for three cases. It is seen that the present prediction for 
m is ~35% larger than the Bartlett et al. prediction. These larger values are the 
result of two combined effects: The first is the larger q c (and q c $) value predicted in 
the present approach. Because these cases are all in the diffusion-limited oxidation 


Table 1: Elemental mass frac- 

tions assumed for the Bartlett et 
al. Avcoat study. 


Element 

Char 

Pyrolysis 

C 

1.0 

0.37 

H 

0.0 

0.13 

O 

0.0 

0.47 

N 

0.0 

0.03 

Si 

0.0 

0.00 
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regime, it can be written from Eq. (3) that 


rh ^ ( 


Qc 


H T -h x 


-)B' 


( 12 ) 


where B' is a constant for this case equal to ~0.4. Although this equation is explicitly applied only in the 
Partially-Coupled approach, it is still approximately true for the Fully-Coupled approach. The second cause 
for the difference between the present results and those of Bartlett et, al. is the influence of treating unequal 
diffusion at the wall. This effect is indicated by the ratio between the Partially- Coupled and Fully-Coupled 
m, which is shown in the last column of the table. The present results are seen to predict a Fully-Coupled 
rh slightly larger than the Partially-Coupled value. Conversely, Bartlett et al. predicts the Fully-Coupled rh 
to be 15 — 20% lower than the Partially-Coupled value. 

Table 2: Comparison of Bartlett et al. results (in parenthesis) with the present Fully-Coupled results (not 

in parenthesis). 


Case 

Uinf 

(m/s) 

Pinf 

(kg/m 3 ) 

Qc, 0 

(W/cm 2 ) 

Qc 

(W/cm 2 ) 

rhxlO 3 

(kg/(m 2 s)) 

T w 

(K) 

™fullv * 

771 partially 

1 

8626 

1.40e-4 

80.7 (66.5) 

70.1 (58.7) 

7.22 (5.34) 

1929 (1990) 

1.02 (0.85) 

2 

7469 

1.87e-4 

61.5 (51.5) 

53.7 (46.4) 

7.42 (5.41) 

1805 (1877) 

1.02 (0.81) 

3 

6099 

2.81e-4 

41.4 (35.4) 

35.6 (33.1) 

7.69 (5.55) 

1633 (1726) 

1.07 (0.76) 


Bartlett et al. values assumed equal to Cm/Ch 


The differences observed in Table 2 are not understood at this time. Possible explanations involve the 
updated thermodynamic and transport properties applied in the present case, or the inaccuracies introduced 
by the bifurcation diffusion approach applied by Bartlett. Also, the present model considers thermo chemcial 
nonequilibrium throughout the flowfield, while Bartlett et al. assumed chemical and thermal equilibrium. 
The nonequilibrium influence was, however, examined by applying LAURA’s free-energy minimization capa- 
bility to Case 1. The computed rh and q c were within 3% of the nonequilibrium values presented in Table 2. 
Thus, the disagreement is not likely due to nonequilibrium effects. 

B. Comparison with Previous Stardust Results 

Olynick et al. 26 and Gupta et al. 2 ' present coupled ablation- flowfield analyses 
for the Stardust vehicle. The study by Olynick et al. coupled a nonequilibrium 
Navier-Stokes flowfield solver to the FIAT material response code. The FIAT code 
applies Approximation #2 to compute ra c , while modeling the in-depth mate- 
rial decomposition to compute rh g . The study by Gupta coupled an equilibrium 
viscous-shock-layer model to a steady-state ablation model. Both of these studies 
treated the pre- flight best estimate trajectory. 

The present results were obtained using the elemental mass fractions for the 
Stardust ablator, PICA, presented by Gupta 27 and listed in Table 3. Steady- state 
ablation is assumed to compute m g and T w . To allow a meaningful comparison 
with Olynick et al., 26 who did not assume steady-state ablation, trajectory points 
were chosen for which Olynick et al.’s rh g and rh c approached the steady-state relationship of Eq. (10), 
therefore indicating near steady- state ablation. The chosen trajectory points are t = 54 and 76 s. Tables 4 
and 5 compare the present stagnation point results with those of Olynick et al. and Gupta. The present 
results apply the Partially-Coupled approach to be consistent with the other studies. For t = 54 s, Table 4 


Table 4: Comparison of Stardust results at the Table 5: Comparison of Stardust results at the 

stagnation point for £ = 54 s. stagnation point for t = 76 s. 



Present 

Olynick et al. 

Gupta 


Present 

Olynick et al. 

Gupta 

9c, o 

920 

1100 

980 

Qc, 0 

206 

205 

200 

Qc 

660 

640 

610 

Qc 

195 

202 

120 

rh 

5.30e-2 

9.00e-2 

3.95e-2 

rh 

3.76e-2 

3.70e-2 

1.95e-2 

T w 

3430 

3400 

3300 

T 

J- -W 

2481 

2550 

2100 


Table 3: Elemental mass frac- 
tions assumed by Gupta for 
PICA. 


Element 

Char 

Pyrolysis 

C 

1.0 

0.60 

H 

0.0 

0.11 

O 

0.0 

0.25 

N 

0.0 

0.04 

Si 

0.0 

0.00 
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shows significant disagreement between the three results. This case is in the sublimation regime, as indicated 
by the large T w values, and is therefore very sensitive to T w and q c . The disagreement in q c ,o suggests that 
the disagreement between the coupled ablation results is caused by differences in flowfield parameters related 
to q c (such as chemical rates or transport properties). For t — 16 s, Table 5 shows good agreement between 
the present results and those of Olynick. This case is in the diffusion-limited oxidation regime, which is the 
regime of interest for the present study. 

V. Analysis with Fixed Ablation Rates 

To evaluate the differences between the Uncoupled, Partially- Coupled, and Fully-Coupled approaches 
defined in Section III, cases with specified m c , rh g , and T w are studied here for each approach. The validity 
of Approximation #1 is assessed by comparing the Ch values predicted by the Uncoupled and Partially- 
Coupled approaches. Recall that Ch is computed from Eq. (1) for the Uncoupled approach, which represents 
Approximation #1, while it is computed directly from Eq. (2) for the Partially- Coupled approach. If the Ch 
values agree between the Uncoupled and Partially-Coupled approaches, then Eqs. (3) and (4) will predict 
the same c w ^ values at the surface (because B' c and B' g will be identical and Cm = Ch )• Thus, if m c was 
to be computed from Eq. (5) instead of held fixed (as will be done in the next section), then the Uncoupled 
and Partially-Coupled approaches would predict the same ra c . However, if the predicted Ch values differ, 
then the predicted m c would be expected to differ accordingly. Generally, a larger Ch translates to a larger 
m c . 

The validity of Approximation #2 is assessed by comparing the c w ,k values predicted by the Partially- 
Coupled and Fully-Coupled approaches. Recall that Approximation #2 is the only difference between these 
two approaches. For a given temperature and pressure, the c w ^ values define the atomic carbon mass 
fraction at the wall (c^c)- The c w ,c values are of interest because they are present in the char equilibrium 
constraint of Eq. (5), which will be applied in the next section to compute m c . Generally, a larger 
in the present analysis will translate to a smaller rh c in the computed-ra c analysis of the next section. It 
will be convenient to compare c w ^ values, as they provide a compact means for interpreting the chemical 
composition at the surface. Typically, larger values of c w ,c are associated with larger values of c w ,o, because 
CO is the dominant species (along with N 2 ) at the surface. The lower the c^o value, the less CO that can 
be formed, and the more c w ,c that is available for forming atomic carbon. 

To examine the influence of the char and pyrolysis elemental composition on the coupled ablation behav- 
ior, the three different compositions listed in Table 6 are considered. The first of these, labeled Avcoat A, 
represents the heritage Apollo Avcoat defined by Bartlett. 25 The second of these, labeled Avcoat B, is the 
same as Avocat A, except with all the H removed. This case will allow the influence of H to be assessed. 
Finally, Avcoat C represents the composition studied in Section IV, which approximates the influence of 
coking and mechanical silica removal. 29 


Table 6: Elemental mass fractions applied in this study for Avcoat. 



Avcoat A 

Avcoat B 

Avcoat C 

Element 

| Char 

Pyrolysis 

| Char 

Pyrolysis 

| Char 

Pyrolysis 

C 

0.49 

0.55 

0.49 

0.60 

1.0 

0.37 

H 

0.0 

0.09 

0.0 

0.0 

0.0 

0.13 

0 

0.27 

0.34 

0.27 

0.38 

0.0 

0.47 

N 

0.0 

0.02 

0.0 

0.02 

0.0 

0.03 

Si 

0.24 

0.0 

0.24 

0.0 

0.0 

0.0 


The present analysis was applied to a 3.6 m radius sphere at a free-stream velocity and density of 
6.0 km/s and 3.0e-4 kg/m 3 , respectively. The specified ablation rates and wall temperatures along the 
body are shown in Fig. 6. So that these specified values would have reasonable profiles and magni- 
tudes, they were obtained assuming steady-state ablation and using the Partially-Coupled approach (for 
Avcoat A). For the laminar case, the resulting stagnation-point surface properties predicted using the Un- 
coupled, Partially-Coupled, and Full-Coupled approaches are listed in Table 7. Results are presented for 
the three Avcoat compositions listed in Table 6. As mentioned previously, the influence of Approximation 
#1 is seen by comparing the Ch values predicted by the Uncoupled and Parti ally- Coupled approaches. 
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For each of the Avcoat compositions, Table 7 shows that 
the Uncoupled Ch is 7 - 15% larger than the Partially- 
Coupled value. This indicates that the Uncoupled ap- 
proach will likely predict, for the computed-m c analy- 
sis presented in the next section, a larger ra c than the 
Partially-Coupled approach. Also mentioned previously, 
the influence of Approximation #2 is seen by comparing 
the c w ,k values predicted by the Partially-Coupled and 
Fully- Coupled approaches. For each of the Avcoat com- 
positions, Table 7 shows that the Parially-Coupled c w ,c 
value is 10% smaller than the Fully- Coupled value, while 
the c Wi o value is roughly 20% smaller. The larger amount 
of oxygen relative to carbon for the Fully-Coupled ap- 
proach results in less atomic carbon at the wall, as seen 
in the last column of Table 7. For the computed -rh c anal- 
ysis, this indicates that the Fully-Coupled approach will 

likely predict a larger rn c than the Partially-Coupled approach. Note that this analysis suggests that Ap- 
proximations #1 and #2 will have offsetting influences in the computed-ra c analysis when comparing the 
Uncoupled and Fully- Coupled m c values. 



Figure 6: Specified ablation rates and wall temperature. 


Table 7: Comparison of Uncoupled, Partially-Coupled, and Fully-Coupled approaches with fixed rh and T w at the stagnation-point. 


Avcoat 

Approach 

Qc 

Cjjx 10 2 

c w ,c 

w,0 

C W ,H 

Cxv,C 

Model 


(W/cm 2 ) 

(kg/(m 2 s)) 





A 

Uncoupled 

38.0 

2.11 

0.171 

0.263 

1.62e-2 

4.1e-17 

A 

Partially- Coupled 

34.6 

1.92 

0.181 

0.264 

1.71e-2 

3.2e-16 

A 

Fully-Couped 

36.4 

1.94 

0.205 

0.322 

1.02e-2 

8.2e-19 

B 

Uncoupled 

38.4 

2.11 

0.180 

0.269 

0.0 

5.4e-14 

B 

Partially- Coupled 

33.4 

1.84 

0.197 

0.272 

0.0 

l.le-11 

B 

Fully-Couped 

34.7 

1.83 

0.207 

0.327 

0.0 

6.9e-19 

C 

Uncoupled 

36.5 

2.11 

0.219 

0.243 

2.26e-2 

9.1e-13 

C 

Partially- Coupled 

34.0 

1.97 

0.229 

0.243 

2.35e-2 

5.7e-13 

C 

Fully-Couped 

35.1 

2.01 

0.277 

0.285 

1.83e-2 

8.2e-13 


The trends mentioned in the previous paragraph for the stagnation-point, are shown to be true down- 
stream of the stagnation point (s = 2 m) in Table 8 for laminar flow. Although only the results for Avcoat A 
are shown, the trends are similar for the other cases, as they were in Table 7. The results for turbulent 
flow are listed in Table 9. The Uncoupled and Partially-Coupled Ch values along the body are compared in 
Fig. 7 for the turbulent and laminar cases. The Ch differences between the Uncoupled and Partially-Coupled 
approaches seen here for the turbulent case are not consistent with the downstream laminar results. Instead 
of the Uncoupled Ch being 7 - 15% larger than the Partially-Coupled value, it is seen to be roughly 15% 
smaller than the Partially-Coupled value. This difference will have a significant influence on the computed 
vn c analysis presented in the next section. This is especially true because the trends in the values 
predicted by the Parti ally- Coupled and Fully-Coupled approaches, and compared along the body in Fig 8, 
are similar to the laminar case. This means that Approximations #1 and #2 will both result in rh c increases 
(in contrast to their offsetting increase and decrease for the laminar case). 


Table 8: Comparison of Uncoupled, Partially-Coupled, and Fully-Coupled approaches with fixed rh and T ^ for the laminar case at s = 2 m. 


Avcoat 

Model 

Approach 

Qc 

(W/cm 2 ) 

ChxIO 2 

(kg/(m 2 s)) 

c w> c 

C'lxj > O 

C\v , H 

c w> c 

A 

Uncoupled 

28.4 

1.56 

0.167 

0.262 

1.58e-2 

l.le-20 

A 

Partially- Coupled 

24.9 

1.38 

0.181 

0.264 

1.71e-2 

8.4e-18 

A 

Fully-Couped 

26.6 

1.40 

0.199 

0.322 

9.74e-3 

1.5e-26 
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Table 9: Comparison of Uncoupled, Partially-Coupled, and Fully-Coupled approaches with fixed to and T w for the turbulent case at s = 2 m. 


Avcoat 

Model 

Approach 

Qc 

(W/cm 2 ) 

Ci/xlO 2 

(kg/(m 2 s)) 

c w ,c 

Civ,0 

Cxu,H 

Cxv,C 

A 

Uncoupled 

35.9 

2.01 

0.198 

0.266 

1.81e-2 

2.3e-13 

A 

Partially- Coupled 

41.5 

2.32 

0.180 

0.264 

1.64e-2 

3.5e-14 

A 

Fully-Couped 

42.3 

2.29 

0.198 

0.309 

l.lle-2 

4.6e-24 



Figure 7: Predicted Cjj values along the body. 



Figure 8: Predicted c u,,k values along the body for the turbulent case. 


The inaccuracy of Approximation #1 in turbulent regions was observed in the previous paragraph. This 
behavior is anticipated to have a significant influence on the computed rn c analysis presented in the next 
section. Therefore, it will be studied in more detail here. The non-ablating (<Z c ,o)> Uncoupled, and Partially- 
Coupled q c values are compared in Fig. 9 for the turbulent Avcoat A case. This figure shows that the 
Partially-Coupled q c becomes nearly equal to the non-ablating value in the downstream turbulent regions. 
This behavior is not followed by the Uncoupled approach, which is the reason for the disagreement in the 
Ch values shown in Fig. 7. Note that q c may be separated into two components, the conductive component: 


Qc,cond 


_ / . dT ve ( 7 

— fcye j + fct 
dz 


dTtr 

dz 


(13) 


and the diffusive component: 

Qc,diff = P ^ ^ (1^) 

i=species 

These components are presented in Fig. 10 for the non-ablating and Partially-Coupled cases. It is seen 
that the conductive component is reduced and the diffusive component is increased with the introduction 
of ablation. The increase of the diffusive component at the stagnation point (therefore, not influenced by 
turbulence) is a result of the different chemistry near the wall for the ablating case. The larger increase of the 
diffusive component downstream of the stagnation point, however, is a result of turbulence. This turbulent 
downstream region was shown in Fig. 4 to depend strongly on the Sc tU rbi and less so on the other details of 
the turbulence model. If Sc tur b is increased above the present value of 1.0, the influence of turbulent diffusion 
decreases. If a large Scturb value is applied to both the non- ablating (which drives the Uncoupled result) 
and Partially-Coupled solutions, the agreement between the Uncoupled and Partially-Coupled q c (and Ch) 
values becomes significantly better. 
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case. 



Figure 10: Comparison of the diffusive and conductive components 
for the turbulent case. 


VI. Analysis with Computed Char Ablation Rates 

The present section engages Eq. (5) to allow m c to be computed as part of the flowfield solution. The 
differences in the ra c predicted using the Uncoupled, Partially-Coupled, and Fully-Coupled approaches are of 
particular interest, and are the main results of this paper. The m g and T w values are fixed to the same values 
as in the previous section. Note that the restriction of fixed T w could be easily removed by applying Eq. 
(9). This was not done here to simplify the interpretation of the results, and because in the diffusion-limited 
regime studied here, T w has a weak influence on ra c . It is assumed that for practical applications a material 
response code will provide rh g , and this value will likely not be sensitive to m c . Thus, holding rh g fixed in 
the present study is justified. 

The differences in <? c , Ch , and c w ^ predicted by the Uncoupled, Parti ally- Coupled, and Fully-Coupled 
approaches were presented in the previous section for fixed m c values (with m g and T w also fixed). In the 
discussion of those results, the results of the present section were correctly anticipated, as will be shown. 

A. Case 1 

The first case to be studied is the same as that considered in Section V, that being a 3.6 m radius sphere 
at a free-stream velocity and density of 6.0 km/s and 3.0e-4 kg/m 3 . The m g and T w values are fixed to 
those shown in Fig. 6. The resulting rn c and q c values are listed in Table 10 for the laminar case at s = 
0.0 m and 2.0 m, and the turbulent case for s = 2.0 m. The m c values computed for Avcoat B are presented 
along the body in Figs. 11 and 12 for the laminar and turbulent cases, respectively. For the laminar case, 
the Fully-Coupled m c is only slightly larger than the Uncoupled result, while it is 10-20% larger than the 
Partially-Coupled result. This peculiar result that the Uncoupled approach is in better agreement (with the 
Fully-Coupled approach) than the Partially-Coupled approach was to be expected from Section V, which 
showed that for the laminar case Approximations #1 and #2 result in offsetting errors. In other words, 
removing Approximation #1 causes the m c difference seen between the Uncoupled and Partially-Coupled 
approaches, which for the laminar case is a 20% decrease in m c . This is a result of Approximation #1 
over-predicting C#, as was shown in Fig. 7. Similarly, removing Approximation #2 causes the m c difference 
seen between the Partially- Coupled and Fully-Coupled approaches, which for the laminar case is about a 
23% increase in rh c . 

For the turbulent results listed in Table 10 and shown in Fig. 12, the comparison between approaches in 
the downstream region is significantly different than the laminar results. As was discussed in Section V, in 
regions of turbulence Approximation #1 under-predicts C#, which is the opposite of the laminar trend. As 
a result, removing Approximation #1 causes an increase in ra c . The influence of Approximation #2 is the 
same for the turbulent case as it is for the laminar case, and therefore removing it also causes an increase 
in rin c . Because removing Approximations #1 and #2 both provide increases in ra c , instead of providing 
offsetting differences as in the laminar case, the difference between the Uncoupled and Fully-Coupled results 
is significant. For the Avcoat B case, the Fully-Coupled m c is 49% larger than the Uncoupled value at 
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s = 2.0 m. 

To examine the influence of the fixed rh g magnitude on the computed ra c , the rn g values applied in 
Table 10 were multiplied by 0.33. The resulting rh c and q c values are listed in Table 11 for the laminar and 
turbulent cases. The trends in the three approaches are similar to those observed previously in Table 10. It is 
therefore concluded that these fundamental trends are not sensitive to m g (although the percent differences 
are). 

Table 10: Comparison of Uncoupled, Partially-Coupled, and Fully-Coupled approaches with fixed rh and for 
the stagnation-point of Case 1. 




s = 0.0, 

Laminar* 

s = 2.0 m, 

Laminar* 

| s = 2.0 m, 

Turbulent 1 

Avcoat 

Model 

Approach 

m c xl0 3 

(kg/(m 2 s)) 

Qc 

(W/cm 2 ) 

rh c xl0 3 

(kg/(m 2 s)) 

Qc 

(W/cm 2 ) 

rh c xl0 3 
(kg/(m 2 s)) 

Qc 

(W/cm 2 ) 


A 

Uncoupled 

5.61 

37.5 

4.21 

27.9 

4.87 

36.9 

A 

Partially-Coupled 

4.87 

34.6 

3.48 

24.9 

6.08 

41.5 

A 

Fully-Couped 

5.53 

35.0 

4.24 

25.3 

6.90 

41.8 

B 

Uncoupled 

5.33 

38.5 

4.01 

28.5 

4.54 

38.0 

B 

Partially-Coupled 

4.27 

34.3 

3.25 

25.5 

5.21 

40.5 

B 

Fully-Couped 

5.64 

33.4 

4.32 

24.7 

6.76 

40.2 

C 

Uncoupled 

3.81 

38.5 

2.80 

28.6 

3.70 

37.3 

C 

Partially-Coupled 

3.46 

35.0 

2.46 

25.2 

4.12 

41.4 

C 

Fully-Couped 

3.68 

36.3 

2.76 

26.2 

4.40 

42.4 


* rhg = 5.48e-3 kg/(m 2 s) 
+ rhg = 3.91e-3 kg/(m 2 s) 
x rhg = 6.33e-3 kg/(m 2 s) 




Figure 11: Predicted laminar rh c values for Case 1 with Avcoat B Figure 12: Predicted turbulent rh c values for Case 1 with Avcoat B 


Table 11: Comparison of Uncoupled, Partially-Coupled, and Fully-Coupled approaches for rh g values equal to 33% 
of those applied in Table 10. 




05 

II 

o 

o 

Laminar* 

s = 2.0 m, 

Laminar* 

| s = 2.0 m, 

Turbulent 

Avcoat 

Approach 

to c x10 3 

Qc 

m c xl0 3 

Qc 

rh c xl0 3 

Qc 

Model 


(kg/(m 2 s)) 

(W/cm 2 ) 

(kg/(m 2 s)) 

(W/cm 2 ) 

(kg/(m 2 s)) 

(W/cm 2 ) 

A 

Uncoupled 

8.52 

38.4 

6.29 

28.5 

8.22 

37.9 

A 

Partially-Coupled 

7.81 

35.6 

5.68 

26.1 

9.22 

41.7 

A 

Fully-Couped 

8.85 

34.8 

6.62 

25.5 

10.4 

41.1 


* rhg - 1.83e-3 kg/(m 2 s) 

* rhg = 1.30e-3 kg/(m 2 s) 
x rhg = 2.11e-3 kg/(m 2 s) 
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B. Case 2 


The second case studied is the same as that considered 
in Section III.D, that being a 3.6 m radius sphere at a 
free-stream velocity and density of 8.0 km/s and 3.0e- 
4 kg/m 3 . Thus, it is the same as Case 1, except with a 
higher velocity. The specified T w and m g are shown in 
Fig. 13, and the resulting rh c and q c values are listed in 
Table 12. The main difference between these results and 
those of Case 1 is that the Partially- Coupled and Fully- 
Coupled results are in closer agreement for Avcoat A and 
C. This is also seen in Figure 14, which shows the rrt c 
distribution along the body for the turbulent Avcoat A 
case. The Avocat B result presented in Figure 15 shows a 
significant difference between the Partially- Coupled and 
Fully- Coupled approaches, similar to that observed for 
Case 1. The difference between Avcoat A and B is that 
Avcoat B contains no hydrogen. It is apparent from the 
present case that the presence of hydrogen reduces the 

disagreement caused by Approximation #2. This is not surprising considering the significantly different 
diffusion characteristics of hydrogen relative to the other heavier species. 

Table 12: Comparison of Uncoupled, Partially-Coupled, and Fully-Coupled approaches with fixed rh and T m for 

the stagnation-point of Case 2. 



Figure 13: Specified m g and T w for Case 2. 




s = 0.0, 

Laminar* 

s = 2.0 m, 

Laminar* 

| s = 2.0 m, 

Turbulent x 

Avcoat 

Approach 

m c xl0 3 

9c 

m c xl0 3 

9c 

m c xl0 3 

9c 

Model 


(kg/(m 2 s)) 

(W/cm 2 ) 

(kg/(m 2 s)) 

(W/cm 2 ) 

(kg/(m 2 s)) 

(W/cm 2 ) 

A 

Uncoupled 

6.90 

84.0 

5.24 

63.3 

6.94 

92.7 

A 

Partially-Coupled 

6.66 

82.3 

4.85 

60.6 

8.99 

106.8 

A 

Fully-Couped 

6.81 

82.7 

5.25 

60.9 

9.02 

106.2 

B 

Uncoupled 

6.49 

85.9 

4.95 

64.4 

6.40 

95.3 

B 

Partially-Coupled 

5.33 

77.9 

4.10 

58.4 

7.21 

100.5 

B 

Fully-Couped 

6.69 

75.9 

5.18 

56.9 

9.19 

99.2 

C 

Uncoupled 

4.87 

86.0 

3.64 

64.9 

5.36 

94.0 

C 

Partially-Coupled 

4.79 

84.6 

3.45 

61.6 

6.21 

108.6 

C 

Fully-Couped 

4.68 

86.7 

3.53 

63.8 

5.95 

108.9 


* m,g = 7.49e-3 kg/(m 2 s) 
+ rh g = 5.46e-3 kg/(m 2 s) 
x mg = 9.39e-3 kg/(m 2 s) 



s (m) 


Figure 14: Predicted turbulent rh c values for Case 2 with Avcoat A. 



Figure 15: Predicted turbulent rh c values for Case 2 with Avcoat B. 
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VII. Concluding Remarks 


A coupled ablation capability was developed for the LAURA Navier-Stokes flowfield code. This capability 
includes the ability to treat a thermochemical nonequilibrium flowfield with turbulence and Stefan-Maxwell 
multicomponent diffusion. A chemical rate model for treating 32 species, including Si species, was compiled. 
The ability to compute m c as part of the flowfield solution was developed for three different approaches. 
The Uncoupled approach, consistent with typically applied design approaches, is simply a post processing 
step to a non-ablating flowfield solution. This approach contains Approximation #1, which represents the 
blowing correction approximation defined in Eq. 1, and Approximation #2, which represents the diffusion 
approximation defined in Eq. 3. The Partially-Coupled approach treats a flowfield with coupled ablation, 
which removes Approximation #1, although it applies Approximation #2. The Fully-Coupled approach 
treats a flowfield with coupled ablation and does not apply either approximation. 

For diffusion-limited oxidation cases relevant to the Orion heatshield composed of Avcoat, the Fully- 
Coupled approach was found to predict rh c values up to 50% larger than the Uncoupled approach. This 
large disagreement was found only for turbulent cases, for which Approximations #1 and #2 both result in 
a lower rn c value. For laminar cases, Approximations #1 and #2 were found to have offsetting effects on 
the predicted rh c . Thus, the Uncoupled and Fully-Coupled results agreed within 10%, while the Partially- 
Coupled value disagreed by as much as 20%. 
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